Spectroscopic Temperature Determination of Degenerate Fermi Gases 
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We suggest a simple method for measuring the temperature of ultra-cold gases made of fermions. 
We show that by using a two-photon Raman probe, it is possible to obtain lineshapes which reveal 
properties of the degenerate sample, notably its temperature T. The proposed method could be 
used with identical fermions in different hyperfine states interacting via s-wave scattering or identical 
fermions in the same hyperfine state via p-wave scattering. We illustrate the applicability of the 
method in realistic conditions for 6 Li prepared in two different hyperfine states. We find that 
temperatures down to 0.05 Tf can be determined by this in-situ method. 

PACS numbers: 



I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sation (BEC) in atomic samples Q have stimulated a 
new wave of theoretical and experimental studies of de- 
generate systems in the dilute and weakly interacting 
regime. Several recent experiments have successfully 
reached the degenerate regime for ultra-cold atomic sam- 
ples of fermions. For example, temperatures of 0.5 Tp 
or less, where Tp is the Fermi temperature of the cold 
fermions, have been reported for 40 K 0, or 6 Li 0. One 
of the goals of reaching even lower temperatures is the 
formation of Cooper pairs of interacting fermions and the 
detection of the Bardeen-Cooper-Schrieffer (BCS) phase 
transition in a trapped gas of fermionic atoms. Such tran- 
sitions in atomic gases will allow the study of phenomena 
usually associated with condensed matter physics, such 
as superconductivity. Recently, "resonance superfluid- 
ity" |4| has been predicted in the strong-coupling regime 
occurring near a Feshbach resonance. This regime is be- 
ing explored experimentally |5j . 

One of the key parameters in studying these systems is 
the accurate value of the temperature T. Various meth- 
ods are being employed to determine T, such as fitting 
the tail of the fermion energy distribution 6] , using mix- 
tures with bosons Q, or using impurities as probes 0. 
Each of these approaches have limitations. For exam- 
ple, considering that the thermal contribution in a Fermi 
gas vanishes with temperature, fitting the tail of the en- 
ergy distribution becomes very challenging. Even in a 
fermion-boson mixture with large attractive interactions, 
such as 40 K- 87 Rb [10, or 6 Li- 7 Li 0, it is questionable 
whether bosons and fermions stay in thermal equilibrium 
when the temperature is lowered to degenerate regime. 
Even if so, the determination of the bosons temperature 
is not easy below T c /2 (T c : BEC critical temperature) 
|10|. Typically, temperature is read by imaging the ve- 
locity distribution of a freely expanded sample. 

In this paper, we describe an in situ method to mea- 
sure the temperature of a sample based on spectroscopic 
measurements. It offers the potential for very accurate 
determination of extremely low temperatures in degener- 
ate Fermi systems as well as other properties (e.g., their 



Fermi temperature Tp). Furthermore, if used in systems 
made of fermions only, it avoids the problem arising from 
the superfluidity of Bose-Einstein condensates when T is 
below the BEC critical temperature T c , which suppresses 
scattering 11]. The approach is based on Doppler-free 
two-color spectroscopy. Such techniques have been suc- 
cessfully employed to create ultra-cold dimers in ultra- 
cold bosonic gases and in BEC |T^ . 



II. METHOD 

The schematic of the two-photon transition is illus- 
trated in Fig. ^ together with the relevant quantities. 
While two atoms approach each other along the molec- 
ular ground state with a given relative kinetic energy e, 
two co-propagating lasers can induce a two-photon tran- 
sition into a bound level of the molecular ground 
state. The first laser of intensity I\ is detuned from the 
excited molecular bound level E^i by 8\ = E^i — v\. Sim- 
ilarly, the second laser of intensity I2 is detuned by an 
amount 62 = -E&2 — {v\ — ^2) from the bound level energy 
E\,2 of the ground molecular state. The rate coefficient 
for the photoassociation (PA) process is given by [l3| 



K(T,Li J i 2 ) = ( ^5^(21 + l)|5 I (e,ii,ia) 



(1) 



where Li = {Ii,Si} represent the parameters of laser i, 
e = h 2 K 2 /2[i = [iv 2 1 2, fi is the reduced mass, and v is 
the relative velocity of the colliding pair. In Eq. (JIJ the 
sum goes over contributing partial waves I, and Si rep- 
resents the scattering matrix element for producing the 
final state E^ from the initial continuum state. Aver- 
aging over relative velocities is implied b y ( ■■•)■ The 
scattering matrix is well approximated by |Rj 
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where 71 is the width of the intermediate bound level 
E b i and 7 S ~ ^ 2 {h/ c)\{E bl \D{R)\e,l)\ 2 is the stimu- 
lated width from the continuum initial state |e, I) to the 
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intermediate state |J3&i). Here, D(R) is the molecular 
dipole transition moment, and in the zero-energy limit, 
7 S oc e 1 / 2+i . Also, 7 = 71 +7 S — li if the laser intensities 
are not too large (this is the regime we are interested in). 
Finally, A± represent a split in the single resonance due 
to the second laser, and is given by 

a± = h& + s 2 ) ± ^{5, - s 2 y + umi , (3) 

where h 2 Sl\ = (2irI 2 /c)\(E bl \D(R)\E b2 )\ 2 defines the 
Rabi frequency f^2- From the analytical form of \Si\ 2 , 
we expect two peaks located at e ~ A±, and a mini- 
mum (| Si | 2 = 0) located between the peaks at e = S 2 
(due to destructive interference between the two scat- 
tering paths). Because of the e 1 ^ 2+ ' Wigner's threshold 
behavior of j s , the peak at A_ will be weaker than the 
one at A + . 

We consider samples of identical fermions of mass m 
in the ultra-cold regime (so that 2// = m). While there 
is no s-wave (I — 0) scattering between two identical 
fcrmionic atoms in the same hyperfine state (the lowest 
contribution is then a p-wave (I — 1) which vanishes as 
e — > 0), s-wave scattering occurs between different hy- 
perfine states. We will concentrate our treatment to the 
latter case, and assume two identical fermions in hyper- 
fine states labeled 1 and 2, respectively. 

To compute the rate coefficient K, we need to aver- 
age over the distribution of relative velocities in the sys- 
tem. The velocity distribution of each ensemble follows 
a Fermi distribution characterized by its temperature T 
and Fermi energy Ep = ksTp. Because of the setups 
in most experiments, we consider the atoms trapped in 
an harmonic potential. For large traps with frequencies 
u)\ , lu 2 , and W3, the density of state for the hyperfine state 
i is given by p(fa) = h 3 k 3 /m 3 uiiuj 2 ui^. 

The Fermi energy Ep of the atoms in a particular hy- 
perfine state is the energy of the highest occupied state 
in the harmonic potential at absolute zero and so is de- 
termined by the total number of atoms sharing the same 
hyperfine state. To have the same number of atoms in 
each hyperfine state provides two benefits: both group 
of atoms have the same Fermi energy, and the s-wave 
scattering rate is maximal. This requirement in not too 
stringent experimentally: if the sample is to be prepared 
from the atoms in a single hyperfine state by transfer- 
ring ~ 50% of the atoms to another hyperfine state, 
then the maximal number difference between the two is 
determined by the degeneracy of the trap potential at 
the Fermi energy. E.g., for isotropic harmonic trap with 
the degeneracy of the m-th energy level ~ m 2 it follows 
AN/N k -iV" 1 / 3 . In what follows we assume that both 
hyperfine states have identical number of atoms so that 
their Fermi energies are identical as well. 

The distribution of the relative momentum re (normal- 



ized to unity) is 
ff 12 (re) = J d 3 fa d 3 k 2 g(fa) g(fa) 5 (fa -fa- 2k) , 

= j d 3 k 2 g(fa) g(k 2 + 2k) , (4) 

where g(fa) = p(ki)fp-D(ki,Ep,T)/Ni. Here, the func- 
tion f FU (fa,E F ,T) = {1 + exp[(££. - E F )/k B T}}- 1 is 
the Fermi-Dirac distribution with E^ — ti 2 k 2 /2m, and 

N{ = J d 3 fap(ki) fpi)(ki 1 Ep, T) is the number of fermions 
in the hyperfine state i. Note that (712 (re) depends on 
Ep and T, but not on the direction of re: the inte- 
gration removes that dependence. In Fig. [21 we illus- 
trate the distribution 312 in term of the relative energy 
e for various temperatures and compare it with the cor- 
responding Maxwell-Boltzmann (MB) distributions. As 
the temperature drops under Tp, the difference between 
the two types of distribution becomes more apparent; as 
the system becomes more degenerate, the gi 2 distribu- 
tion spreads over a larger range of energies (when com- 
pared to MB distributions). As T approaches zero, 1712 
is entirely contained between e = and 2Ep. As op- 
posed to the single particle distribution g(E) (see inset), 
which is contained between E — and Ep at zero-T and 
decreases sharply to zero at Ep, g\ 2 goes to zero more 
gradually. In fact, for the case of free particles, one can 
show that 312(e) is a triangle with values 1 at e = and 
at e = 2Ep, while g(E) is the standard step function 
between E = and Ep. 

We want to use 2-photon scattering to probe the rela- 
tive velocity (or energy) distribution and infer the tem- 
perature of the system. At ultra-cold temperatures, 
Eq.JIJ becomes K = (-Kh/2m) J d 3 K 3i2(re)|S , | 2 /re, which 
can be rewritten in terms of e as 

A 2 rco 

K(T,L 1 ,L 2 ) = —J deg 12 (e)\S (e,L 1 ,L 2 )\ 2 . (5) 

Besides T and Ep, K depends on the detunings 8\ and 5 2 
and the Rabi frequencies 1 and ^2- Because of the im- 
portant variation of |<So| 2 with respect to S 2 , lineshapes 
obtained as a function of 5 2 contain precious information 
about the degenerate Fermi gas. As T is reduced, the 
distribution 312 becomes more sharply defined for e be- 
tween and 2Ep, and by scanning S 2 for nearby values, 
the large variations in the integrand of Eq. JSJ will probe 

312- 

We computed lineshapes for realistic systems. The 
Fermi energies Ep reached in todays experiments are in 
the range of few /iK. In addition to this experimental 
constraint, one needs to be detuned far enough from reso- 
nance to not populate the level v\ , with typical detunings 
Si - 150 MHz (7.2 mK) [3. The lifetime of most levels 
vi is about 10-20 nsec, giving 71 - 16 - 32 MHz (0.76- 
0.38 mK) Finally, one needs intensities large enough to 
drive the transition, but not too large to heat the sys- 
tem (through broadening). Typical values for high lying 
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states v\ at 500 GHz below the asymptote are fii = 100 
kHz (4.8 LiK) and Q 2 ~ 2-3 MHz (96-144 ^K). Under 
these conditions, we have Si ^> hfl 2 . With large Si, A + 
will also be large, and the only way to get a peak of IS"! 2 
with e ~ Ep will come from A_: hence, we must have 
62 ~ Ep as well. Then, Si 3> hfl 2 ^> S 2 and we will have 
peaks at A_ ~ S 2 — h?Vt^/Si and A + ~ S\ + ft, 2 f2 2 /(5i, as 
well as a minimum at e = 5 2 . 

III. RESULTS AND DISCUSSION 

We obtained lineshapes for various temperatures by 
varying S 2 in j£j . We selected parameters consistent with 
the case of 6 Li (although similar results can be obtained 
for other species): Ep = 1 /aK, Si = 7 mK, 71 = 0.76 
mK and hQ 2 = 100 /iK. The remaining parameters lead 
to a multiplicative constant in K that can be factored out 
for all lineshapes. We show representative lineshapes nor- 
malized to unity (in arbitrary units) in Fig. [31 as T is low- 
ered below Tp, the peak becomes narrower and sharper. 
Many other features can be identified on those lineshapes. 
As 5 2 grows, K reaches a minimum value K m i n , increases 
rapidly to its maximum value K max , and then decreases 
more slowly to a "background" value Kb g . In addition, in- 
flection points on both sides of the peak can be identified. 
All these features can be extracted from the lineshapes, 
and could be used to extract temperature. For example, 
one could take the ratio [K max — Kb g ]/[K m i n — Kb g ] of a 
given lineshape, in principle removing the multiplicative 
constant in K. However, this ratio varies very little for 
T/Tp ~ 0.5 and below. Besides the effect of a back- 
ground signal/noise would also affect/reduce the sensi- 
tivity of the temperature determination. One could also 
map the location (i.e. the value of 5 2 ) of these features 
as a function of T, but again, their variations are small 
below say T/T F ~ 0.1. 

Alternative features, such as ratios of the derivatives 
at the inflection points (on each side of the main peak) 
could be used. In this case, the background would not 
play any role, nor the multiplicative constant in K. A 
simple measure incorporating both criteria is the ratio 
of the curvature of the signal at the peak value and its 
height as defined from the minimum value of K (see inset 
in Fig. i.e. 

C^- K K ( V El. (6) 

While finding the second derivative of a lineshape with 
respect to 5 2 may be challanging in general, at the local 
peak maximum this should not be a problem as the first 
derivative there is zero. 

Because of the rapid increase of K from its minimum 
value to its peak value, only a small range of S 2 needs 
to be scanned (as compared to the ratio of the deriva- 
tives at the inflection points): for Ep ~ 1 ^K, scans 
with kHz precision are required. In Fig. 01 we show C 
as a function of T /Tp: it decreases monotonically with 



T, and even at T/Tp ~ 0.05, its variation is noticeable. 
So, by measuring three features on the lineshape, namely 
its minimum, maximum, and the curvature at the max- 
imum, one can evaluate C and determine the tempera- 
ture of the degenerate Fermi sample by reading it from 
Fig. Although the exact values of C depend on other 
parameters (71, 5±, and il 2 ), its general shape follows the 
example shown in Fig. 0] For illustration purposes, we 
also include two other curves (corresponding to other sets 
of Si and f^)- The curve with the largest variations at 
low T /Tp is found when we have <5i 3> hfl 2 3> S 2 with a 
factor of 100 between each parameter. 

The lineshapes shown in Fig.|3are not very sensitive to 
small variations in 71, <5i or £l 2 : variations of few % give 
very similar curves. Of these parameters, 71 cannot be 
varied experimentally, except by changing the intermedi- 
ate molecular level En- However, for most atoms with 
Fermi species considered in todays experiments, i.e., al- 
kali metals, 71 varies little with En (e.g., for 6 Li2 in the 
triplet excited state, 71 varies from ~ 11.7 MHz down 
to 3.2 MHz (or 0.56 mK to 0.16 mK). However, systems 
with even smaller 71 (about 1-10 /iK) would offer the pos- 
sibility of measuring the velocity distribution gi 2 directly 
by scanning S 2 . In fact, \S\ 2 then becomes proportional 
to the delta-function S(e — S 2 ). The extremely small vari- 
ations in gi 2 as T — * could possibly be detected. An- 
other possibility to enhance the sensitivity of this type 
of measurement could be realized using Feshbach reso- 
nances: a sharp increase in the integrand of Eq.Q at 
e ~ e ros , where e rcs is the resonant relative energy, would 
emphasize the corresponding region of the velocity distri- 
bution. By positioning the resonance appropriately (e.g., 
via external magnetic fields), strong contribution of re- 
gions of gi2 where small variations occur as T — * could 
be obtained, hence more precise information could be 
extracted from the lineshapes. Notice that p-wave Fes- 
hbach resonances 0] could be used to enhance p-w&ve 
scattering in a single-fermion population. 



IV. CONCLUSION 

We have shown that the lineshapes obtained using a 
two-photon Raman probe can yield very precise deter- 
mination of the temperature in degenerate Fermi gases. 
One of the key advantage of this approach, beside its 
non-destructive nature, is the extremely low tempera- 
tures measurable: in the example considered here, as low 
as ~ 0.05TV. Although we have illustrated this two-color 
scheme for identical fermions in two different internal 
states with equal populations, this method can be made 
much more flexible. For example, if the populations are 
very different, the gi 2 distribution may exhibit sharper 
structures that could be used advantageously. The en- 
hancement due to Feshbach resonances could also help 
probing chosen regions of gi 2l and extend the method 
to samples of identical fermion in a single internal state, 
by using p-wave Feshbach resonances. Finally, when the 
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interaction between the fermions is attractive such as 
for 6 Li, Cooper pairs may form. Photoassociation of a 
Cooper pair should demand less energy than the pho- 
toassociation of the two fermions produced by breaking 
either one or two Cooper pairs 16]. These processes 
should contribute additional features on top of the sim- 



ple lineshapes, their locations and magnitudes probing 
the energy gap and the number of Cooper pairs. 
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FIG. 2: Relative velocity distribution for various temperatures: gi2(K)d 3 K = 2n(2jj,/h 2 ) 3 ^ 2 gi2(e)y' r ede (here, <?i2(e)\/e is 
normalized to unity). The effect of degeneracy becomes apparent when gi2 (solid lines) are compared to the corresponding 
Maxwell-Boltzmann distributions (dashed lines). The inset shows the corresponding single-particle distribution g(E)\fE. 
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FIG. 3: Lineshapes as the functions of 82 for temperatures T/Tf = 0.1, 0.5, and 0.9, with 71 = 3 mK, Si — 0.56 mK, and 
0,2 — 100 n~K. The definitions of K", K max , K min , and (5 2 nax are illustrated in the inset. 
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FIG. 4: Ratio C from Eq.@ as a function of T/Tf- The sensitivity of C on other parameters is also illustrated by two other 
cases (dashed lines). 



